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We present a microscopic model for a strongly repulsive electron gas on a 2D square lattice. We 
suggest that nearest neighbor Coulomb repulsion stabilizes a state in which electrons undergo a 
"somersault" in their internal spin-space (spin-flux). When this spin- 1/2 antiferromagnetic (AFM) 
insulator is doped, the charge carriers nucleate mobile, charged, bosonic vortex solitons accompa- 
nied by unoccupied states deep inside the Mott-Hubbard charge-transfer gap. This model provides 
a unified microscopic basis for (i) non-Fermi-liquid transport properties, (ii) mid-infrared optical ab- 
sorption, (iii) destruction of AFM long range order with doping, (iv) angled resolved spectroscopy 
(ARPES), and (v) d-wave preformed charged carrier pairs. We use the Configuration Interaction 
(CI) method to study the quantum translational and rotational properties of such pairs. The CI 
method systematically describes fiuctuation and quantum tunneling corrections to the Hartree-Fock 
approximation and recaptures essential features of the (Bethe ansatz) exact solution of the Hubbard 
model in ID. For a single hole in the 2D AFM plane, we find a precursor to spin-charge separation. 
The CI ground state consists of a bound vortex-antivortex pair, one vortex carrying the charge and 
the other one carrying the spin of the doping hole. 
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I. INTRODUCTION 

In 1986 Bednorz and MuUeifl discovered that the 
perovskite (BaLa)2Cu04 exhibits high-temperature su- 
perconductivity, with a critical temperature of up to 
30K. Soon after, La2-:ESrj:Cu04 and YBa2Cu307_:r were 
found to have superconducting critical temperatures of 
35K and 95K respectively. Since then, many such com- 
pounds were found, including the Tl and Hg series. The 
current record Tc of 135K (165K under pressure) is found 
in the HgBa2Ca2Cu308 system. A typical phase diagram 
is shown in Fig. |l} The undoped parent is an insulator 
with long range antiferromagnetic order. Extremely low 
doping ( a; « 0.02 charge carriers per site) leads to a 
complete destruction of the long-range AFM order, and 
a transition to an unusual non-Fermi-liquid metal. This 
unusual metal becomes superconducting, with the transi- 
tion temperature Tc strongly dependent on the doping x. 
The maximum Tc is reached for dopings around x — 0.15. 
For higher dopings the critical temperature decreases to 
zero, and in the overdoped region a crossover towards a 
(non-superconducting) Fermi-liquid takes place. 

The effective two-dimensional Hamiltonian we use 
to describe the electrons residing in the 0(2po-)- 
Cu(3da;2_j^2) orbitals of the isolated Cu02 plane is the 
generalized one-band Hubbard Hamiltonian: 

n = -tY^ {t,jc\^Cja + h.c}j + V.jn.rij (1) 

where cj^ creates an electron (in the orbital centered) at 
site i with spin tr, tij is the hopping amplitude from site j 

to site i on the square lattice, Ui = ^la'^i'y total 
number of electrons at site i, and Vij is the Coulomb re- 
pulsion between electrons at sites i and j. The dominant 
terms are the nearest-neighbor hopping tij — to and the 
on-site Coulomb repulsion Va = U/2. If only these two 
terms are considered, and we shift the chemical potential 
by U, this reduces to the widely studied Hubbard model. 
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FIG. 1. Schematic phase diagram as a function of doping 
of high-temperature cuprate compounds. 



In the undoped parent compound phase, there is one 
electron per orbital, and in the absence of interactions 
one would expect these compounds to be metallic, with 
a half- filled conduction band. Instead, from Fig. we 
see that they are insulators with long-range AFM or- 
der. This is a strongly-correlated electron system, with a 
large on-site repulsion term { U > t), which causes elec- 
trons to become localized one per each orbital in order 
to avoid energetically expensive double occupancy. This 
obviously leads to insulating behavior, while the antifer- 
romagnetism is simply, a perturbational effect from the 
(small) hopping term.cl However, as the plane is doped 
with holes, some of the previously filled orbitals are emp- 
tied, and electrons in nearby orbitals can freely hop into 
them, allowing for charge conduction. This is also seen 
from the phase diagram in Fig.^, which shows that a very 
small amount of doping x ~ 0.02 completely destroys the 
LR AFM insulator phase, and the system becomes metal- 
lic (or superconductor, at lower temperatures). However, 
this metal has very unusual non- Fermi liquid properties. 
Further doping leads to a crossover to a more conven- 
tional (but non-superconducting) Fermi liquid metal in 
the extremely overdoped region. 

Understanding the non-Fermi liquid metal above the 
superconducting state is the central issue in the cuprate 
physics. Some of the most striking evidence of non- 
Fermi-liquid behavior is provided by Angle Resolved 
Photo-Emission Spectroscopy (ARPES), which clearly 
shows the absence of quasi-particle peaks in the normal 
state. Equally compelling evidence is provided by resis- 
tivity measurements p which reveal a scattering rate in- 
versely proportional to the temperature r ~ l/r, extend- 
ing over a range of up to TOOK.! In an ordinary Fermi- 
liquid, electron-electron scattering gives a dependence 
of the relaxation rate, related to the quadratic energy de- 
pendence of the quasiparticle lifetime i ^ |e — e_Fp. This 
is a hallmark of a Fermi liquid B In fact, the canonical 
behavior is indeed observed in the extremalv overdoped 
region, already identified as a Fermi-liquid.l3 But its ab- 
sence in the intermediate doping region of the unusual 
metal, combined with the absence of quasiparticle peaks 
in the ARPES data clearly show that this unusual metal 
is not a Fermi-liquid. Then, a natural question arises. 
If the charge carriers of the unusual metal are not the 
quasiparticle-like charge carriers of a Fermi-liquid, what 
is their nature? 

Some clues are provided by experiments. Hall mea- 
surements in the underdoped regime tell us that the 
charge carriers are positively charged aad their density 
equals the doping, i.e. the hole density.Q (Deviations of 
both the sign and the density scaling are observed for 
some compounds in the overdoped regions). Optical mea- 
surements reveal the appearance, with doping, of a low- 
frequency Drude tail, which suggests the existence of free, 
or very mobile charge carriers J3 These optical measure- 
ments confirm both the anomalous 1 /T scattering rate of 
these charge carriers, as well as the fact that their density 
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equals the doping. Together, these measurements suggest 
that each hole introduced in the Cu02 layer evolves into 
a mobile, positively charged, carrier (but not a quasipar- 
ticle). However, both types of measurements show other 
anomalous behavior as well. The Hall coefficient has a 
strong, 1/T temperature dependence, which is very puz- 
zling, since neither the density nor the charge of the carri- 
ers is temperature dependent, and neither is the Lorentz 
force. The optical measurements reveal that the Drude 
tail is accompanied by the appearance of a broacL mid- 
infrared band deep inside the Mott-Hubbard gap.B This 
mid-infrared band develops with doping, and signifies the 
existence of localized levels inside the Mott-Hubbard gap. 
Unlike the Drude tail which collapses below Tc to a 5{Ld) 
function, this mid-infrared response persists unchanged 
both above and below Tc. Magnetic measurements indi- 
cate that although the LR AFM is completely suppressed 
above x « 0.02, strong short-range AFM correlations 
persist, with a correlation length roughly equal to the. 
average distance between the holes (charge carriers). □ 
Neutron scattering indicates the appearance of incom- 
mensurate peaks in the magnetic structure factor, with 
a shift from the AFM (7r/a,7r/a) vector varying linearly 
with doping for 0.02 < x < 0.12 and then saturating. 
The four incommensurate peaks are arranged diagonally 
for X < 0.05 and then rotate by 45° for x > 0.05.0 The 
doping dependence of all these features clearly indicates 
that magnetism is crucial to the entire phenomenology. 

More puzzling behavior appears in the superconduct- 
ing state. Flux quantization clearly proves that pairing 
does take place, and the unit of charge in the super- 
conducting condensate is 2e. ARPES and phase sensi- 
tive measurements have shown that the superconducting 
gap has d-wave symmetry, and this is believed to mirror 
the symmetry of the "Cooper-pair" wave-function to in- 
ternal rotations. Penetration depth measurements show 
that the density of superfluid in the limit T — > scales 
with doping. This means that the "Cooper-pairs" must, 
in fact, be formed from the positive charge carriers of 
the unusual metal, not from electrons, as in conventional 
BCS theory (the density of electrons is 1 — x, not x). This 
underscores the need to identify the charge carrier of the 
unusual metal, but leads to another question, namely, 
how does strong pairing (leading to high superconduct- 
ing temperatures) occur in a system dominated by strong 
Coulomb repulsion? 

In this paper we describe a microscopic model which 
offers simple and compelling answers to the above ques- 
tions, as well as to other puzzling features described 
above. Unlike other approaches which assume that the 
fundamental quantum degrees of freedom of the many- 
electron system are conventional and that the resulting 
phenomenology is an "emergent law of nature" arising 
from the complexity of the system, we propose that there 
is a hidden fundamental law of Nature. This fundamen- 
tal law of nature expresses a novel dynamical degree of 
freedom, namely that an electron can perform a "somer- 
sault" in its internal space of Euler angles (when consid- 



ered as a rigid body of non-zero volume) as the electron 
traverses a closed loop in external coordinate space. The 
result is a new quantum number for the many-electron 
system (the eigenvalues of the spin- flux), which we pro- 
pose is as fundamental as the existence of the spin-i 
itself. We argue that the electronic somersault (spin- 
flux) is dynamically generated through electromagnetic 
interactions and in particular, the off-diagonal part of 
the Coulomb repulsion between electrons. In Section || 
we introduce the spin-flux Hamiltonian, which we pro- 
pose as the appropriate Hami lton ian describing the iso- 
lated Cu02 plane. In Section HI we apply Hartree-Fock 
Approximation (HFA) and the configuration interaction 
(CI) method to the spin-flux Hamiltonian and we iden- 
tify both the mobile bosonic charge carrier as well as the 
nature and symmetry of the strong pairing interaction 
between such charge carriers. We provide a physical jus- 
tification for the accuracy of the CI approximation and 
we explicitly demonstrate this by comparison with the 
exact solution of the ID Hubbard model. t3 Finally, in 
Section IV we discuss the comparison of our model and 



its results with the experimental findings and draw some 
final conclusions. 



II. THE SPIN-FLUX MODEL 

The neglect of the dynamical consequences of longer 
range Coulomb interaction {Vij = O.ifi ^ j), in the 
generalized Hubbard model of Eq. (|l|), is based on the 
assumption of uniform charge distribution and on the 
Fermi-liquid theory notion of screening of the effective 
electron-electron interaction. However, Fermi liquid the- 
ory fails to explain many of the crucial features of the 
high-Tc cuprates. In our description, we include the 
nearest-neighbor Coulomb repulsion, which we assume is 
on the energy scale of t. This has important dynam- 
ical consequences in our model and cannot simply be 
absorbed into the Madelung constant. In particular, it 
leads to the generation of spin-flux, an entirely new type 
of broken symmetry in the many-electron system, which 
we show leads naturally to bosonic charge carriers in the 
form of meron-vortices, non-Fermi-liquid behavior and a 
strong attractive pairing force between holes in the AFM 
background. 

The concept of spin-flux is closely related to the ex- 
istence of the spin— i particles in nature. The spin 
of a physical electron may be regarded as arising fro™, 
the quantization of a classical, symmetric spinning toptil 
whose kinematical properties are described by a set of 
three independent Euler angles. These Euler angles con- 
stitute a continuous manifold, the group manifold of 
SO (3). Unlike the manifold S2 (the surface of a unit 
sphere) describing the orientation of a classical magnetic 
moment, the group manifold of SO (3) is not topologically 
simply connected. According to the axioms of quantum 
mechanics any physical wave-function must be every- 
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where continuous and difFerentiable. For this to be sat- 
isfied on a simply connected manifold, the wave-function 
must be single valued. Consequently the 0(3) nonlinear 
sigma model describes integer spins. The doubly con- 
nected group manifold of SO (3), however, can accom- 
modate two- valued wave-functions which are everywhere 
continuous and differcntiable. This leads to half-integer 
spins. In order to accommodate spin— i, the 0(3) nonlin- 
ear sigma model must be supplemented with a magnetic 
monopole which is placed at the center of the sphere S2.t3 
The charge of this monopole corresponds to a quantum 
of the third Euler angle in the parameterization of the 
SO (3) group manifold. 
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FIG. 2. Conventional (a) vs. spin-flux (b) trajectory. 
In the latter case, the spin rotates by 2n (somersault) as it 
encircles the path in real space. 

Spin-flux arises when the electron executes a topolog- 
ically nontrivial path within its internal space of Euler 
angles while it traverses a closed loop in the external 
coordinate space. The internal path is one that cannot 
continuously be deformed to zero and corresponds to a 2tt 
rotation (somersault) in the space of Euler angles. An- 
other depiction of this process (see FigJ|) is seen by con- 
sidering distinct points within the spinning top and fol- 
lowing their trajectories as the electron executes a closed 
path in external coordinate space. In a non-spin-flux cir- 
cuit, the two trajectories are unlinked, whereas in a spin- 
flux circuit (involving a somersault) the two trajectories 
are linked. We propose that this novel possibility repre- 
sents a hidden but fundamental Law of Nature which has 
not been considered in conventional treatments of many- 
body theory. Spin-flux corresponds to a fundamentally 
new quantum number in a many-electron system and re- 
quires the extension of the conventional many-electron 
Hilbert space. We suggest that this simple addition (at a 



fundamental level) to the dynamical degrees of freedom 
of interacting electrons an a two-dimensional lattice leads 
to a unifled, microscopic explanation of a large variety of 
experiments on the cuprates. 

In order to describe the above physics from our starting 
Hamiltonian, 

Ti^-to^2 (cl^Cja + h.c^+U^^nii^n^i + V^n^nj (2) 

we introduce bilinear combination of electron operators 
Kj = 4a'^ai3(^30' M = 0,1,2,3, for i ^ j (summation 
over multiple indexes is assumed). Here fi" is the 2x2 
identity matrix and a = (a^,a^,a^) are the usual Pauli 
spin matrices. The notation means that the sites 

i and j are nearest neighbors. The quantum expecta- 
tion value ( ) of the A^^ operators are associated with 
charge-currents (/i — 0) and spin-currents (/x = 1,2,3). 
Non- vanishing charge currents lead to appearance of elec- 
tromagnetic fields, which break the time-reversal sym- 
metry of the Hamiltonian. Experimentally, this does not 
occur in the cuprates. In the following, we adopt the 
ansatz that there is no charge current in the ground state 
A^^ — but circulating spin-currents may arise and take 
the form Af_j- = ^P-iAijna,a = 1,2,3, where \Aij\ = A 
for all i and j, and n is a unit vector. These spin-currents 
provide a transition state to the uniform spin-fiux mean- 
field that we use in this paper. In principle, non-uniform 
states of spin- flux may arise, in which | Ay | has a nontriv- 
ial deppH-dence on i and j. One such case was discussed 
earlier ,Ej in which skyrmion textures in the AFM back- 
ground carry quantized units of spin-flux. In this case 
Ay is a dynamical variable. However, for the purpose of 
this paper, we consider only a uniform, static, mean-fleld 
conflguration of the spin-flux. 

Using the Pauli spin-matrix identity, |o'^^(cr^,^,)* = 
Saa'Sfsp', it is possible to rewrite the nearest-neighbor 
electron-electron interaction terms as n^rij = 2ni — 
^A^(A^)+. We neglect fluctuations in the spin- 
currents, and use the mean-field factorization to replace 
A^(A^)+ (A^)(A^)++A^;(A^;)*-(A^;)(A^)*. Thus, 
the quartic nearest-neighbor Coulomb interaction term is 
reduced to a quadratic term that is added to the hopping 
term leading to the effective Hamiltonian: 

^ ^ E {'^I^L^jP + '^■c-) + E "»T"«i- (3) 

Here, T^jj = ((5^^ + iA^ti ■ aa(i)l\J\ + A^ are spin- 
dependent SU (2) hopping matrix elements defined by 
the mean-field theory, and t = ioVT+A^. In deriv- 
ing (^ we have dropped constant terms which simply 
change the zero of energy as well as terms proportional 
to ^ rii which simply change the chemical potential. It 

was shown previouslyEUli that the ground state energy 
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of the Hamiltonian of Eq.(||) depends on the SU(2) ma- 
trices T*-' only through the plaquette matrix product 
^12^232^342.41 = exp(m • ct$). Here, $ is the spin-flux 
which passes through each plaquette and 2<I> is the angle 
through which the internal coordinate system of the elec- 
tron rotates as it encircles the plaquette. Since the elec- 
tron spinor wave- function is two- valued, there are only 
two possible choices for $. If $ = we can set T^jj = 6ij 
and the Hamiltonian (^) describes conventional ordered 
magnetic states of the Hubbard model. The other pos- 
sibility is that a spin-flux $ = tt penetrates each pla- 
quette, leading to T^^T'^3rp34rp4i ^ _i j.^^^ means 

that the one-electron wave-functions are antisymmetric 
around each of the plaquettes, i.e. that as an electron 
encircles a plaquette, its wave-function in the internal 
spin space of Euler angles rotates by 27r in response to 
strong interactions with the other electrons. In effect, the 
electron performs an internal "soaaersault" as it traverses 
a closed path in the Cu02 plane.E3 This spin-flux phase is 
accompanied by a AFM local moment background (with 
reduced magnitude relative to the AFM phase of the 
conventional Hubbard model). In the spin-flux phase, 
the kinetic energy term in exhibits broken symmetry 
as though a spin-orbit interaction has been added. In 
the presence of charge carriers this mean-field is unstable 
to the proliferation of topological fluctuations (magnetic 
solitons) which eventually destroy AFM long range order. 
In this sense, the analysis which we present below goes 
beyond simple mean field theory. The quantum dynam- 
ics of these magnetic solitons described by the Configu- 
ration Interaction (CI) method, corresponds to tunneling 
effects not contained in the Hartree-Fock approximation. 
For simplicity, throughout this paper we assume that the 
mean-field spin-flux parameters T''^ are given by the sim- 
plest possible choice T^^ ^ -l^T'^^ = T^-^ ^ T^^ = 1 
(for more details, see Ref. 14). In order to go beyond 
a mean-fleld description of the spin-flux, these matrices 
may also be treated as dynamical variables. In this paper, 
we go beyond mean-field theory in describing the antifer- 
romagnetic degrees of freedom but restrict ourselves to a 
mean-field model of the spin-fiux. 

The mean-field ground-state of the undoped spin-flux 
model is an AFM Mott insulator. It is interesting to 
note that the quasi-particle disjjifirsion relation obtained 
in the presence of the spin-fluxEj accurately recaptures 
the dispersion as measured through angle-resolved photo- 
emission spectroscopy (ARPES) in a compound such as 
Sr2Cu02Cl2li3 (see Fig. ||). There is a a large peak cen- 
tered at (7r/2,7r/2) with an isotropic dispersion relation 
around it, observed on both the (0, 0) to (tt, tt) and (0, tt) 
to (tt, 0) lines. The spin-flux model in HFA exhibits an- 
other smaller peak at (0, 7r/2) which has been observed 
in more recent experimental data£3 The quasi-particle 
dispersion relation of the conventional Hubbard model 
(ri2 = = = T^i = 1) has a large peak at 
(7r/2,7r/2) on the (0,0) to (tt, tt) line (see Fig. |), but it 
is perfectly flat on the (0, tt) to (tt, 0) line (which is part 



of the large nested Fermi surface of the conventional 2D 
Hubbard model). Also, it has a large crossing from the 
upper to the lower band-edge on the (0, 0) to (0, tt) line. 
Both this dispersion relation and the very similar one of 
the t — J model (see Ref. 15) are in contradiction to 
ARPES data. 
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FIG. 3. A comparison between the experimentally de- 
termined E{k) quasi-particle dispersion relation, from angle 
resolved photoemission studies (ARPES), for the insulating 
Sr2Cu02Cl2 and the HP AFM spin- flux model dispersion rela- 
tion (full line) and the HP AFM conventional Hubbard model 
dispersion relation (dashed line). Three directions in fc— space 
are shown: (0,0) to (tt.tt), (7r,0) to (0,0) and (7r,0) to (0,7r). 
While the peak on the (0, 0) to (tt, n) is equally well described 
in both models, the mean-field spin-flux model gives a much 
better agreement for the (7r,0) to (0,0) and (7r,0) to (0, tt) 
directions. The fitting corresponds to U = 2.01 eV, t — 0.29 
eV for the spin-fiux phase, and U = 1.98 eV, t — 0.21 eV 
in the conventional phase. The experimental results are the 
ARPES results of Ref. 15 (circles) and Ref. 16 (triangles). 



III. DOPING INDUCED MERON- VORTEX 
SOLITONS 

A. The Static Hartree-Fock Approximation 

The HP results for the undoped AFM ground state 
of the spin-flux Hamiltonian are in good agreement with 
experimentally measured dispersion (see Fig. ^ . The az- 
imuthal symmetry of the dispersion relations about the 
Fermi points plays a key role in determining the symme- 
tries of doping induced magnetic configurations. This is 
more straightforward to see in a simpler, continuum ver- 
sion of the model, obtained by letting the lattice constant 
a ^ (see Ref. 17,18). Since the dispersion relation near 
the Fermi point q = (^, ^) is isotropic, the dependence 

on k — K — q ^ — iVf of the continuum HP equations 
is such that it preserves rotational invariance. As a re- 
sult, the 2D HP equation reduces trivially to a ID radial 
equation, with a structure very similar to that of the ID 
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difFerentiai HF equation obtained for the ID Hubbard 
model.EZItJ Once this radial ID solution is found, the 2D 
configuration is simply generated through a 2n rotation 
about an axis perpendicular to the 2D plane. As a result, 
there is a close analogy between solutions obtained for the 
ID Hubbard modeJ^-aiid Jor the 2D spin-flux model, in all 
our investigations .ll2lllZril2l 
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FIG. 4. Self-consistent spin distributions of a 10x10 lattice 
with a meron- vortex in the spin-flux phase. The core of the 
meron is localized in the center of a plaquette, in the spin-flux 
phase (in the conventional phase, the core of the meron- vortex 
is localized at a site). This excitation has a topological wind- 
ing number 1, since the spins on either sublattice rotate by 
2-K on any curve surrounding the core. The magnitude of the 
staggered magnetic moments is slightly diminished near the 
vortex core but is equal to that of the undoped AFM back- 
ground far from the core. The hole is localized on in the 
vortex core. 

In the corresponding discrete model, the relevant 
doped configuration is the meron- vortex (see Fig.^). The 
doping hole is trapped in the core of a magnetic vortex, 
which indeed has azimuthal symmetry. The b 0.1.1 nd level 
on which the hole is trapped can be showr£3o to be 
split from the top of the valence band and drawn deep 
inside the Mott-Hubbard gap. As a result, the meron- 
vortex is a charged boson. This can be inferred by direct 
inspection of Fig.^, which shows a configuration with to- 
tal spin zero and a positive charge trapped in its core. 
An argument based on the electronic structure, identical 
to the one offered for the charged-bpaMiic domain-walls 
of the polyacetylene, also holds .liZHla'tSEj The parallel to 
the quasi-one dimensional ID polyacetylene is again a 
reflection of azimuthal symmetry which reduces the 2D 
continuum model to a ID radial equation. Clearly, the 
isotropic dispersion about the Fermi points is responsi- 
ble for the appearance of bosonic charge carriers. They 
are very unlike quasi-particle charge carriers, which carry 



both spin and charge together. Thus, one would expect 
a metal with such bosonic charge carriers to have intrin- 
sically different properties from those of a Fermi-liquid 
metal. 
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FIG. 5. Self-consistent spin distribution for a tightly 
bound meron-antimeron pair. The meron (M) and the an- 
timeron (A) are localized on neighboring sites. The total 
winding number of the pair is zero. The magnetic AFM order 
is disturbed only on the region where the vortices are local- 
ized. The attraction between holes is of topological nature 
and on long length scale is stronger than unscreened Coulomb 
repulsion between charges. The doping charge is mostly lo- 
calized on the two plaquettes containing the meron and an- 
timeron cores. The two holes localized in the vortex cores are 
responsible for the fact that the meron-antimeron pair does 
not collapse. 

A look at the spin-configuration in Fig.^ also shows 
that this cannot be realized if cyclic boundary conditions 
(CBC) are imposed. With them, the self-consistent so- 
lution found in the presence of one hole is a very differ- 
ent configuration called the spin-bag, which is a rather 
immobile quasiparticle-Jike, configuration (carrying both 
charge +e and spin-i).L30 However, for more than one 
hole added to the AFM plane, the HF ground-state of the 
spin-fiux Hamiltonian always shows meron-vortices cre- 
ated through doping, even with CBCEJ The simple rea- 
son for this is that while a single isolated meron-vortex 
is incompatible with CBC, meron-antimeron pairs are 
not (see Fig. |^ for a typical meron-antimeron configura- 
tion). In fact, the nucleation of merons and antimerons 
in pairs also solves another problem, related to the en- 
ergetic cost of creating a meron. It is straightforward 
to prove that the energy of a single memn configuration 
increases logarithmically with its size.Ej Beyond a cer- 
tain (large) separation it is energetically favorable for the 
meron-antimeron pair to collapse into a pair of charged 
spin bags. In the conventional Hubbard model (with no 
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spin-flux) spin bags are favored at all separations. In the 
spin-flux phase, the farther the meron is from the an- 
timeron, the more spins in between are rotated by the 
vortices, and the excitation energy increases. As a re- 
sult, an isolated pair tends to be as closely bound as 
possible. If the vortices were uncharged, at low temper- 
atures they would annihilate each other. However, the 
holes localized in the vortex cores lead to a very strong 
short-range Coulomb repulsion which prevents the pair 
annihilation, thus making the pair stable. It is worth 
noting that even in the complete absence of screening, at 
long distances the 1 jr Coulomb repulsion would be over- 
come by the ln(r) attraction between vortices, leading to 
a stable bound pair. 

To conclude, we see that even at the static HF level, 
charge carriers in the spin-flux phase exhibit bosonic na- 
ture, and a strong pairing attraction to other charge 
carriers. This attraction is of magnetic origin, arising 
from exchange energy lost by spins which are no longer 
perfectly AFM aligned. The importance of the bosonic 
meron-vortex excitations becomes even more apparent 
when we consider fluctuation and tunneling corrections 
to the HF approximation. These correspond to transla- 
tional motion of the charged vortices. It turns out that 
charged meron-vortices have an effective mass compara- 
ble to that of the band electron. As seen from the CI 
method, they are much more mobile than the very heavy 
spin-bags. 

B. Fluctuations and Tunneling: The Configuration 
Interaction Method 

We now discuss how to improve the mean-field de- 
scription above. Given a complete basis of states {| </>!)} 
spanning the iV-body Hilbert space, the exact TV-body 
ground state wave- function can be written as j^*) — 
X^jO!*!*^*)- "^^^ coefficients {ai} are found from solving 
the Schrodinger equation = which reduces 

to a linear system of equations Ti-ijaj = Ej^j ^ijO^j^ 
for all i. Here, Hij = {(j)i\H\4ij) and Otj = {4>i\4>j) are the 
matrix elements of the Hamiltonian in the complete basis 
chosen and the overlapping matrix for the basis states, 
respectively. 

In general, the number of basis states of the iV-body 
problem increases exponentially with N, and the prob- 
lem becomes untractable for values of N which are rather 
small. To deal with large values of N, one is forced to 
truncate the complete basis set, and only retain a smaller 
subset of states which are most likely to contribute sub- 
stantially in the decomposition of the ground-state wave- 
function. The choice of this smaller subset is the crucial 
issue, since it determines the quality of the approxima- 
tion. In what follows, we describe a particular subset 
which allows us to recapture certain key features of the 
exact Bethe ansatz solution of the ID Hubbard model. 

Consider the set of Slater determinants generated 



from the Hartree-Fock solution. Its states are \^hf) 
and all possible particle-hole excitations {a\ap\^HF}}, 
{ajjajj/OpOp' I'I'/zf)} etc. Obviously, if all possible com- 
binations of occupied and empty orbitals are considered, 
the set thus spanned is a complete basis of the A^-body 
space. 

Let us now order the states in this HF basis set accord- 
ing to their energies {(j)\'H\(j)) . If we are interested in the 
ground-state and the low-lying excitations of the system, 
we only need to keep the low-energy states of the HF basis 
set. This procedure is, in fact, very well known for Hamil- 
tonians with a non- degenerate HF ground-state. If only 
the HF ground-state and the states with one particle-hole 
excitations are kept, this leads to the Random Phase Ap- 
proximation (RPA). Besides a better approximation for 
the ground-state than simple HFA, the RPA enables us 
to find collective excitations and the particle-hole contin- 
uum. 
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FIG. 6. The lowest energy dispersion band Ej{k) (in 
units of t) as a function of the total momentum k of the 
meron-antimeron pair. The momentum units are n/a and 
U/t = 5. For convenience, the reference energy is taken to be 
the static HF energy of the self-consistent meron-antimeron 
pair. Quantum hopping and rotation lowers the overall energy 
of the pair by 1.75t. 

On the other hand, the doped HF ground-states of 
the spin-fiux Hamiltonian are degenerate. For instance, 
the meron-antimeron pair shown in Fig.^ happens to be 
centered at site (10,10). It is obvious that configurations 
which have the meron-antimeron pair centered about any 
other site will have the exact same HF energy (cyclic 
boundary conditions are imposed). Also, if the center for 
the pair is fixed at a site, there are four distinct possi- 
ble arrangements of the meron and antimeron about that 
site, obtained by rotating the meron-antimeron axis by 
90°. Thus, for a plaquette of size N^xNy, the HF ground 
state is AN^Ny degenerate. Clearly, the minimal choice 
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for the subset of states used to search for the ground-state 
of the doped system must include all these degenerate HF 
ground-state wave-functions. This choice is the essence 
of the configuration interaction (CI) approximation.Ell 
One expects these HF states to mix with equal weight 
jofip in the decomposition of the ground-state and low- 
energy states l^*) = " c'il'Pi)- Specifically, we de- 
note \'^e{n,m)) to be the HF ground-state wave- 
function describing a configuration centered at site (n, m) 
of the lattice, and with the meron-antimeron axis at an 
angle 9 from the x-axis. Then, 1 < n < Nx, 1 < m < Ny 
and e = 45°, 135°, 225° or 315°. On general symmetry 
grounds one expects the ground-state and the low-lying 
energy states to have the general form 



l*j(fc)> 



E 



The cyclic boundary conditions limit the vector fc to a 
subset of equally spaced values inside the first Brillouin 
zone and J must be an integer. 
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FIG. 7. The rotational symmetry of the meron-antimeron 
wave-function as a function of the total momentum carried 
by the pair (measured in units of n/a). The outside region 
(containing the absolute minima points (7r,7r)) has d-wave 
symmetry (J = 2), while the core region about the (0,0) 
point has s-wave symmetry {J — 0). The intermediary area 
is a mix of s+d wave symmetry. 

Eq.(||) shows that the ground-state and low- lying en- 
ergy states found within the CI method have transla- 
tional and rotational symmetry. This procedure over- 
comes the most glaring shortcoming of the mean-field 
theory (the broken translational and rotational invari- 
ance). Clearly, the CI wave- functions describe quantum 
dynamics of the charge carriers. The pair is no longer 
pinned at one site, as in HFA, but moves over the entire 
lattice. The mobility of the pair and its preferred internal 
angular momentum can be obtained from the dispersion 
relation Ej{k) = j{k)\Ti\'ii j{k)) . The lowest energy 
band of the meroruantimeron pair obtained for U/t = 5 
is shown in Fig.^.t2l We find that the CI ground-state 



corresponds to pairs of total momentum (tt/o, n/a), and 
which have d-wave symmetry, J = 2. In fact, J varies 
throughout the Brillouin zone as shown in Figj^, from 
pure d-wave around the (±7r/a, ±7r/a) points to pure s- 
wave around the (0,0) point, where a second local min- 
imum exists. The large width (« At) of the dispersion 
band clearly proves that the meron-antimeron pair is a 
very mobile excitation. The existence of mobile charged 
bosonic merons and antimerons thus provides a micro- 
scopic basis for the non- Fermi liquid "parent" metal from 
which superconductivity emerges. 



HF — Charged domain-wall 
HF - Spin-polaron 
Bedie ansacz result 




FIG. 8. Excitation energy, in units of t, for a mobile 
charged domain wall (squares, dashed line) and a mobile 
charged spin polaron (squares, full line), as obtained from 
the CI approach. The exact excitation energy given by 
the Bethe-ansatz method is shown by diamonds. The do- 
main-wall CI energy is in excellent agreement with the ex- 
act BA results overt the entire U/t range, while the spin po- 
laron CI energy is significantly different. For comparison, we 
also show the excitation energies for the domain wall (cir- 
cles, dashed line) and the spin-polaron (circles, full line) as 
obtained within HFA, proving again that the extra kinetic en- 
ergy gained by the moving domain wall is of order t for most 
U/t values. In contrast, the extra kinetic energy gained by 
the spin polaron is of order t^ /U — > as U/t increases, so in 
the large U/t limit there is almost no difference between the 
HF and CI results for the charged spin polaron. We conclude 
that the charged domain-wall is the relevant excitation for all 
values of U/t. 

As an indication of the validity and accuracy of the CI 
method, we briefly review thejesults of such an analy- 
sis of the ID Hubbard model.E3 Here, the analog of the 
a charged meron-vortex is the charged domain wall soli- 
ton which facilitates a tt flip from one AFM ground-state 
to the other AFM ground-state (with all spins flipped). 
The charged domain-wall traps the hole on a pair of lev- 
els that are localized deep,iiiside the Mott-Hubbard gap, 
and is a charged bosonO'ESM The analog of the spin- 
bag is the spin-polaron, which disturbs the AFM order 
only locally, trapping the hole in its small FM core. As a 
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result, the spin-polaron carries both charge and spin— i. 
In Fig.|| we show the HF (circles) and CI (squares) ener- 
gies of both the spin-polaron (full line) and the charged 
domain- wall (dashed line). For the spin-polaron, the dif- 
ference between the HF and CI energy varies as /U . 
In analogy with the 2D case, this is the expected behav- 
ior since the spin— ^ carried by the spin-polaron localizes 
it on one magnetic sublattice. On the other hand, the 
domain-wall lowers its energy by an amount of the or- 
der t between the HF and the CI values. This clearly 
shows that the charged domain- wall is a very mobile ob- 
jects, and that its dynamics must be properly described 
in order to get a realistic picture. When this is done, 
it is apparent that the charged domain-wall is the low 
energy excitation over the entire U jt range, and in fact 
its energy is in excellent agreement with the exact Bethe 
Ansatz prediction (shown as diamonds). 

Finally, we consider the stability of the ground-state 
and the low-lying energy states to inclusion of more states 
in the subset used to generate them. The next meaningful 
enlargement of this subset is to add to it all states with 
one particle-hole excitation, obtained from all degener- 
ate HF ground-state wave- functions. While this increases 
the size of the subset substantially, it does not lead to a 
lower ground-state energy. This is because the system is 
gapped. The lowest energy particle-hole excitation is ob- 
tained when an electron from the top of the valence band 
is excited on one of the bound levels, localized in the cores 
of the meron and the antimeron. These levels are roughly 
one quarter of the way into the gap (for J7/t = 5) , so the 
states with a particle hole excitation have an energy of 
order U S more than the HF ground-states. Consequently 
they do not mix into the ground-state. They only influ- 
ence the higher energy states. The low-lying excitations 
depicted in Fig. ^ are charge excitations, associated with 
motion of the charged meron-antimeron pair. These are 
not affected by addition of excited states in the varia- 
tional subset, except at rather high energy. But there 
also exist low-lying spin-wave excitations, which are de- 
scribed by adding particle-hole excited states (RPA) in 
the variational subset. While it is very likely that the 
nucleation of merons and antimerons in the AFM back- 
ground will alter the dispersion of the spin-waves, it is 
interesting to note that low-lying charge and low-lying 
spin excitations have very different origin. This is re- 
lated to the spin-charge separation tendency exhibited 
by this system. 

If a single hole is added to the 2D AFM plane and the 
cyclic boundary conditions are imposed, the HF approx- 
imation leads to a spin-bag mean-field solution for the 
spin-flux phase. We analyze the translational properties 
of the spin-bag using the CI method. As it turns out,t£l 
is that the spin-bag is a very immobile object. The width 
of its dispersion band is of the order of t^/J7, i.e. much 
smaller than the ~ t bandwidth of the meron dispersion. 
The reason is that the spin-bag carries both charge and 
spin-i. The spin-^ only allows the spin-bag to live on 
one of the magnetic sublattices. (If it were to move to 



the other sublattice, its spin should flip, and such pro- 
cesses are forbidden) . In order to move, a spin-bag must 
tunnel two sites to a second-nearest neighbor site, and 
this is a t'^/U process. On the other hand, the charged 
meron carries no spin, so there are no restrictions for its 
motion. It can move to a nearest neighbor site, leading to 
a t- hopping process. This suggests that when dynamics 
is properly taken into consideration within CI, a highly 
mobile vortex-antivortex pair sharing a single hole is en- 
ergetically favored to the immobile spin-bag. We have 
verified this hypothesis, and demonstrated that an upper 
bound to the energy of a singly-charged vortex-antivortex 
pairjs, indeed, much lower than the CI energy of spin- 
bagJl^ The charged spin-bag solution is thus unstable to 
dissociation into a highly mobile vortex-antivortex pair, 
which shares the hole. When a second hole is added, 
a second single-charged vortex-antivortex pair is nucle- 
ated. However, it becomes energetically more convenient 
for the two holes to become bound to the same vortex- 
antivortex pair, leading to the appearance of a meron- 
antimeron pair which carries both charges (pre-formed D- 
wave pair). The remaining uncharged vortex-antivortex 
pair is unstable to collapse, at low temperatures. 

IV. DISCUSSION AND COMPARISON WITH 
EXPERIMENTS 

The crucial distinguishing feature of our model is the 
concept of spin-flux, the dynamical possibility of an elec- 
tron undergoing a somersault as it traverses a closed loop. 
In the undoped parent compound, this leads to dispersion 
relations with isotropic symmetry about the Fermi point, 
in excellent agreement with those measured experimen- 
tally through ARPES. This symmetry of the dispersion 
relations leads to real-space configurations which have 
the same type of symmetry. In the spin-flux model, the 
holes doped into the AFM plane nucleate magnetic vor- 
tices and become trapped in their cores, leading to the 
appearance of mobile, bosonic, charge carriers. 

In contrast, both the conventional Hubbard model and 
its asymptotic limit, the t-J model, exhibit a very large, 
nested Fermi surface (at the mean-field level) in the un- 
doped parent compound. This Fermi surface has quasi- 
linear ID character, since there is no fc-dependence along 
the Fermi surface. As a result, configurations stabi- 
lized by doping exhibit the same quasi-lD character in 
real space. It has been suggested they take the form of 
charged stripes.cS We note finally that vortex-like config- 
urations are unstable in the conventional Hubbard model, 
while stripe-like configurations are generally unstable to 
the formation of a quantum liquid of merons in the spin- 
flux model. One notable exception is the commensurate 
case X = 0.125, when merons and antimerons crystallize 
along ID lines, leading ta the stripe-like configuration 
observed experimentally,E3 provided that ajSmall (3%) 
hopping anisotropy is included in the model.Efl 
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The mobile bosonic charged meron-vortices created 
through doping provide a microscopic basis for a non- 
Fermi hquid behavior. They also exhibit a very strong 
pairing attraction, of magnetic origin. This pairing, 
which arises in a purely repulsive electron system, leads 
to appearance of pre-formed "Cooper-pairs" of d-wave 
symmetry. This agrees with the experimental findings 
of d-wave superconductivity and the scaling of the su- 
perfluid density with doping. The pre-formed pairs may 
also be related to the observation of a jpseudo-gap on the 
underdoped side of the phase diagram.On3 

Many other features of our model are in agreement 
with experimentally observed properties of the cuprate 
superconductors. Nucleation of magnetic vortices with 
doping explains a variety of magnetic properties, starting 
with complete destruction of the long-range AFM order 
for very low doping concentration. As we can see from 
Fig. ^, a tightly-bound meron-antimeron pairs disturbs 
the long-range AFM ordering of about 100 spins of the 
lattice. For very low dopings, these pairs are far from 
each other, and there are many spins on the plane whose 
orientations are not affected by any pair. Thus, most of 
the spins maintain the long range AFM order. However, 
as the doping increases over about 2%, the areas occu- 
pied by each meron-antimeron pairs start to overlap with 
those occupied by the neighboring pairs. At this doping 
the orientation of all the spins on the Cu02 planes is af- 
fected by at least one pair of vortices, and therefore the 
LRO is lost. The local ordering, however, is still AFM. 
This picture explains the extremely low doping necessary 
for the disappearance of LR AFM order, as well as the 
fact that the spin correlation length is basically equal to 
the average distance between holes (voilices) and does 
not depend strongly on the temperature.l3 Each hole car- 
ries its vortex with it, and the spins in each vortex are 
correlated with each other. The correlation length is thus 
roughly equal to the average inter- vortex (inter-hole) dis- 
tance. The nucleation of magnetic vortices quantitatively 
explains the split of the (tt, tt) AFM Bragg peak into the 
four incommensurate peaks whose positions shift with 
dopingjj as. observed in LaCuO and, more recently, in 
YBaCuO.L^ The form factor of any given vortex already 
gives rise to an apparent splitting of the neutron scatter- 
ing peak. As demonstrated in Ref. 14, even at the mean- 
field level we recapture the neutron scattering data using 
the HF distribution of meron-vorticesr, A more detailed 
investigation suggests that saturationE3 of the peak split- 
ting for X > 0.12 may be related to expansion of tlaacore 
radius of the individual vortices at higher doping.Es 

Optical behavior of the cuprates is also explained nat- 
urally using our model. Two features develop in the 
optical absorption spectra with doping: a broad mid- 
infrared temperature-independent absorption band, and 
a st|Cpngly temperature-dependent low-frequency Drude 
taiLQ In our model the broad mid-infrared band is re- 
lated to excitation of electrons from the valence band 
onto the empty levels bound in the vortex cores ,0 which 
are localized approximately one quarter of the way inside 



the Mott-Hubbard gap. The number of localized levels 
scales with the number of vortices, and inter-vortex and 
spin-wave interactions lead to their broadening into the 
observed band. This mechanism is similar to the one 
leading to a broad mid-infrared absorption band in poly- 
acetylene with dopingo (In the continuum limit, the 
meron vortex in fact creates a pair of mid-gap states 
in the Mott-Hubbard gapEJ). The polyacetylene band 
is due to electronic excitations inside the cores of the 
polyacetylene domain- wall solitons ,Ej which are the topo- 
logical analogues of meron-vortices .EZrIia Another strong 
argument in favor of this interpretatipn is provided by 
photoinduced absorption experiments E3 If the undoped 
parent compounds are illuminated with intense visible 
light, they develop absorption bands that resemble the 
mid-infrared bands of the doped compounds. Similar be- 
havior is observed in polyacetylene, and is attributed to 
the nucleation of solitons by photoexcited electron-hole 
pairs.EiJ 

The second component of the optical spectrum is the 
Drude tail. It results from the response of the freely mov- 
ing charged vortices to the external electric field. The 
strong temperature dependence of this tail is determined 
by the scattering mechanism for merons (due to inter- 
actions with other merons and spin- waves). This inter- 
pretation is also supported by the fact that the super- 
conducting transition leaves the mid-infrared absorption 
band unchanged. Merons with internal electronic struc- 
ture are still present on the planes but pair condensation 
leads to a collapse of the Drude tail into a 5{lo) response. 

Finally, our model provides some understanding of the 
cross-over towards the Fermi-liquid metal in the over- 
doped phase. For large dopings [6 > 0.30 — 0.40) the av- 
erage inter-vortex spacing becomes extremely small and 
the very cores of the merons start to overlap. In this limit 
the Mott-Hubbard gap is completely filled in by the dis- 
crete levels, and the spin-flux state becomes-energetically 
unstable relative to a normal Fermi liquidJlil 

We note, finally, that the spin-fiux Hamiltonian has es- 
sentially no free or adjustable parameters. The choice of 
U/t is fixed by the experimentally measured size of the 
Mott-Hubbard charge transfer gap at zero doping. All 
other experimental features such as (i) the position and 
nature of the mid-infrared optical absorption band, (ii) 
the ARPES data and (iii) the position of the magnetic 
neutron scattering satellite peaks as a function of doping, 
are quantitatively described by the same choice. More de- 
tailed comparisons with specific experiments may require 
the incorporation of specific (smaller energy scale) inter- 
actions which are not included in this simplest version of 
the spin-flux Hamiltonian. 
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